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ABSTRACT 


This paper proposes and analyses a basic deterministic mathematical 
model to investigate Simulation for controlling malaria Diseases 
Transmission dynamics. The model has nine non-linear differential 
equations which describe the control of malaria with two state variables for 
mosquitoes populations and five state variables for humans population and 
to introduce the new SPITR model for the transmission dynamics of malaria 
with four time dependent control measures in Ethiopia Insecticide treated 
bed nets (ITNS), Treatments, Indoor Residual Spray (IRs) and Intermittent 
preventive treatment of malaria in pregnancy (IPTP). The models are 
analyzed qualitatively to determine criteria for control of a malaria 
transmission dynamics, and are used to calculate the basic reproduction/^. 
The equilibria of malaria models are determined. In addition to having a 
disease-free equilibrium, which is globally asymptotically stable when the 
R 0 < 1, the basic malaria model manifest one's possession of (a quality of) 
the phenomenon of backward bifurcation where a stable disease-free 
equilibrium co-exists (at the same time) with a stable endemic equilibrium 
for a certain range of associated reproduction number less than one. The 
results also designing the effects of some model parameters, the infection 
rate and biting rate. The numerical analysis and numerical simulation 
results of the model suggested that the most effective strategies for 
controlling or eradicating the spread of malaria were suggest using 
insecticide treated bed nets, indoor residual spraying, prompt effective 
diagnosis and treatment of infected individuals with vaccination is more 
effective for children. 

KEYWORDS: Malaria, Basic reproduction number, Backward bifurcation 
analysis and Vaccination 


INTRODUCTION 

Malaria is an infectious disease and is life threatening for 
human beings worldwide. Parasite is an organism that 
lives on or inside a human body from which it gets its food. 
Malaria is caused due to a parasite called Plasmodium. 
Plasmodium parasite is transmitted into human body 
when an infected female anopheles mosquito makes bites. 
Plasmodium parasites making the human liver as their 
home multiply their population and start infecting red 
blood cells of the human. A variety of plasmodium 
parasites exist. Mainly four types of plasmodium cause 
malaria disease among the human viz., falciparum, vivax, 
ovale and plasmodium malaria [2], [22], [34] and [35]. 

The plasmodium parasite is injected into the human 
bloodstream in the form or stage or life cycle known as 
sporozoite. The parasites go through a complex life cycle 
inside the hosting human body and they live at various 
stages both in liver and red blood cells. From time to time 
the parasites pass through various stages of their life cycle 
and during which numerous human red blood cells are 
destroyed. 

The effect of malaria disease varies with the infecting 
variety species of Plasmodium and also with prior health 
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and immune status of the individual. Typically malaria 
disease causes fever and chills together with headaches, 
vomiting and diarrhea. It may also cause long-term 
anemia, liver damage and neurological damage. The most 
dangerous falciparum parasite scan causes cerebral 
malaria. It causes frequently a fatal condition involving 
damaging the brain and central nervous system. The 
survived people from the cerebral malaria may too 
experience brain damage. 

Malaria is a parasitic vector or mosquito borne disease. 
The malaria disease is found among particular people or in 
a particular region in many parts of the world [3]. 
Annually about 3 to 5 hundred million malaria cases are 
being identified worldwide and of which 1 to 3 million 
people end their lives with deaths [4]. 

About 40% of the world population lives in malaria 
endemic regions. The parasites and mosquitoes causing 
malaria have been developing their resistance powers 
against the drugs and chemicals that are being used to kill 
or destroy them. As a result, the incidents of malaria 
disease have been increasing in the past few decades. 
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In 2018, 49 malaria-endemic countries reported less than 
10 000 cases of the disease, up from 40 countries in 2010. 
Twenty-seven countries had less than 100 cases - a strong 
indicator that elimination is within reach - up from 17 
countries in 2010. In 2018, two countries were awarded 
the official WHO certification of malaria elimination: 
Paraguay and Uzbekistan. Two more countries Algeria and 
Argentina were certified malaria-free in 2019. The 
certification is granted when a country is able to prove, 
beyond a reasonable doubt, that the chain of indigenous 
malaria transmission has been interrupted for at least 
three consecutive years. Globally, a total of 38 countries 
and territories have achieved this milestone [9]. 

In the recent time, significant resources and control 
programs have been made available worldwide. The aim is 
to reduce malaria infected cases and prevalence or gain 
upper hand over the disease. Different strategies and 
programs with varying effectiveness and efficiency are 
being adopted to control malaria disease. Comparative 
knowledge of these existing programs is necessary to 
design and organize any new and useful and cost effective 
procedure to control malaria epidemic [5]. 

Efforts to reduce malaria transmission have led to the 
development of efficient vector control interventions, 
particularly insecticide treated nets (ITNs), indoor 
residual spraying (IRS), and larvicide [7-8]. The ITNs 
include conventional nets treated with a WHO 
recommended insecticide and long-lasting insecticidal 
nets. Note that larva is an immature form of an insect and 
larvicide is a chemical used to kill larvae. These 
interventions are used in malaria endemic countries 
especially those in sub-Saharan Africa and have led to 
reduction in malaria morbidity and mortality 
substantially. However, malaria epidemic continues to 
claim hundreds of thousands of lives every year, thus 
necessitating a continued control effort to fight against the 
disease [9 - 10]. 

One year ago, Gilimbeta was sitting on a wooden bench at 
a rural health facility on the outskirts of Lilongwe, Malawi, 
with a mixture of excitement and nervousness. That day, 
her daughter, Lusitana, would be the first child in the 
world to receive the RTS malaria vaccine through routine 
immunization as part of a landmark pilot in Ghana, Kenya 
and Malawi. One year later, a total of 275,000 children like 
Lusitana have received their first dose of the RTS malaria 
vaccine through the pilot implementation programme 
[ 10 ]. 

The RTS malaria vaccine is the first and, to date, the only 
vaccine that can significantly reduce malaria in children - 
the group at highest risk of dying from malaria. The 
vaccine is complementary to other proven measures to 
fight malaria. Having grown up in Malawi, a country where 
suffering from malaria is very common, Gilimbeta knows 
all too well the toll the disease takes on children, families 
and communities. 

Since 2000, the scale up of proven interventions, such as 
insecticide-treated bed nets, has driven down malaria 
illness and death, but progress has slowed, particularly in 
Africa. Efforts are underway to increase the use of the 
tools we have, but even with scale up of existing tools 


there is a critical need for new tools to further drive down 
malaria illness and deaths the RTS malaria vaccine may be 
one of those tools. The malaria vaccine can reduce malaria 
in children, including cases of severe malaria, related 
hospital admissions and the need for blood transfusions. 
The current pilot programme will help guide and optimize 
the vaccine's future use. 

Malaria and immunization services critical during COVID- 
19 prevention is especially important at this time. We 
know from experience that a disruption in essential health 
services can seriously increase illness and death from 
common diseases, including malaria and other vaccine 
preventable diseases, particularly among the most 
vulnerable. As the COVID-19 pandemic evolves, 
immunization against vaccine preventable diseases, and 
continued use of malaria prevention measures are critical 
for reducing illness and relieving pressure on fragile 
health systems. 

Vaccines are an essential part of protecting children from 
deadly diseases and should continue to be delivered. 
Immunizations can be delivered in ways that minimize any 
health risks for health workers and communities, to help 
ensure that vaccines reach everyone, everywhere [10-11]. 

Malaria has been considered as a global issue. 
Epidemiologists together with other scientists invest their 
efforts to understand the dynamics of malaria and to 
control transmission of the disease. From interactions 
with these scientists, mathematicians have developed 
tools called mathematical models. These models are used 
significantly and effectively for giving an insight into (i) 
the interaction between the host and vector population (ii) 
the dynamics of malaria disease (iii) control mechanisms 
of malaria transmission and (iv) effectivness of eradication 
techniques. 

Mathematical models are particularly helpful as they 
consider and include the relative effects of various 
sociological, biological and environmental factors on the 
spread of disease. The models have played a very 
important role in the development of malaria 
epidemiology. Analysis of mathematical models is 
important because they help in understanding the present 
and future spreads of malaria so that suitable control 
techniques can be adopted. 

The SIR is a simplest mathematical model and has three 
classes or compartments viz., Susceptible, Infected and 
Recovered. Here we considered SIR model and modified it 
by adding two more compartments viz., protected and 
treatment classes. The inclusion of the protected and 
treatment classes to SIR model and extending it to SPITR 
model is a technique used to control the spread of malaria 
disease. Thus, we presented the SPITR model that 
describes the dynamics and controlling mechanism of 
malaria disease. The effect of controlling technique in the 
spread of malaria is analyzed [12]. This new model is an 
extension or modification of the existing mathematical 
models used to deal with malaria epidemic viz., SIR, SEIR, 
SPEITR, SEIS, SEIIR [12 - 15], [18] and [22]. The 
simulation studies of the model with variable values of 
sensitive parameter of the spread of malaria are 
performed and the results are incorporated. The necessary 
conclusions have been drawn. 
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1. The SPITR model 

Malaria disease occurs due to the bites of female 
anopheles mosquito. That is, the prevention and control 
rates of malaria disease will increase with the mortality or 
death rate of the female Anopheles mosquito. Hence, the 
need of inventing intervention strategies which will help 
to increase mortality of the female anopheles mosquitoes 
in addition to the natural mortality is arisen. In the recent 
times many intervention strategies are adopted 
worldwide. These intervention strategies increase the 
mortality rate of female anopheles mosquitoes and help 
directly to control the birth and spread of malaria. The 
important intervention strategies include (i) indoor 
residual spraying (IRS) and (ii) insecticide treated bed net 
(ITN). These intervention strategies play a significant role 
in the controlling mechanism of malaria disease [1]. 

1.1. Classification of human compartments 

The present mathematical model is constructed with an 
objective of incorporating the effect of intervention 
strategies on the spread of malaria disease. Thus, we 
considered the protected and treatment class, added to 
SIR model and extended to SPITR model. In the SPITR 
model all the humans are divided into five classes or 
compartments and they exhibit the following properties: 

Susceptible human class (S h ) contains humans those do 
not have malaria disease but are likely to be bitten by 
infected female anopheles mosquitoes causing malaria. 
The people in this class do not take any protective 
measurements against malaria disease. In general, all the 
common people are included in this class. Humans of the 
susceptible class have the potential to remain either in the 
same class or to enter into protected class or infected 
class. People in this class are given awareness about 
malaria disease and advised to be protected using IRS and 
ITN mechanisms. The people who follow advises and 
adopt protected mechanisms will enter into protected 
class. The remaining people will remain in this class as 
long as they do not be bitten by infected female anopheles 
mosquitoes and enter into infected class. 

Protected human class (P)contains humans those do not 
have malaria disease. Also this class people are assumed 
that they will never get the disease. The people in this 
class do take some protective measurements like IRS and 
ITN and they become fully protected from malaria disease. 
These people keep themselves safe from malaria. The 
humans of this class are drawn from both susceptible and 
recovered classes. 

Infected human class (I h ) contains humans those are 
already infected and got malaria disease. People come into 
infected class from susceptible class. People from infected 
class have got the potential to enter into either treatment 
class or recovered class. 

Treatment human class (T) contains the humans those 
are already infected and got malaria disease. The people in 
this class are being given medical care or clinical 
treatment. People into the treatment class come from only 


infected class. After completion of the treatment 
successfully, the people from treatment class will go to 
recovered class. 

Recovered human class (R) contains the people who 
recover from the malaria disease andreturn to normal 
status of health. People come into the recovered class from 
both infected and treatment classes. People in this class 
are given awareness about malaria disease and advised to 
be protected using IRS and ITN mechanisms. The people 
who follow advises and adopt protected mechanisms will 
enter into protected class and the remaining people will go 
to susceptible class. 

Further, we assume that the sick people do not travel far 
away or they go from place to place. Thus, we do not 
include the immigrations of malaria infected people in our 
model. 

1.2. Classification of mosquito compartments 

Anopheles is a mosquito of a genus which includes the 
species that transmits the malaria parasite to humans. 
Anopheles male mosquitoes do not carry plasmodium 
parasite carrying capacity and they do not cause malaria in 
humans. Hence, anopheles male mosquitoes are not 
included in the model. Infected female anopheles 
mosquitoes through bits cause malaria in humans. We 
have thus considered and included only female anopheles 
mosquito population in the present model. Further, note 
that only female anopheles mosquitoes bite humans for 
blood meals. In the present model all the female anopheles 
mosquitoes are divided into two classes or compartments 
and they exhibit the following properties: 

Susceptible mosquito class (S m ) contains the female 
anopheles mosquitoes. The mosquitoes in this class are 
not infected and their bites do not lead to malaria disease 
in humans. The population size of these mosquitoes grows 
with the natural birth rate ofpand diminishes with the 
natural death rate of co. The susceptible female anopheles 
mosquitoes get infected on feeding themselves with the 
blood of malaria - infected humans and then they leave 
susceptible mosquito class and enter infected mosquito 
class. Susceptible mosquitoes bite the infected humans 
with the rate of Izl. Further, it is assumed that the 
susceptible mosquitoes get infected with a rate of 
6 hm w hen they bite infected humans. 

Infected mosquito class (/ m )containsthe infected female 
anopheles mosquitoes. The mosquitoes in this class are 
infected and their bites lead to malaria disease in humans. 
The population size of this class grows with a transfer rate 
proportional to a m from the susceptible mosquito class. 
The infected mosquitoes diminish with the natural death 
rate of co. The infected female anopheles mosquitoes 
transfer infection to susceptible humans when the former 
feed them with the blood of the latter. Infected mosquitoes 
bite the susceptible humans with the rate of 0. Further, it 
is assumed that the infected mosquitoes transfer infection 
with a rate of 6 mh w hen they bite susceptible humans. 
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YW HUMAN POPULATION 



MOSUITO 

POPULATION 

Figure 1 Flow chart of SPITR malaria model 

The parameters in Table 1 and the state variables in Table 2 are used in Figure 1 to formulate the malaria model with 
control. 


Parameter 

Description 

¥ 

Natural birth rate of humans 

P 

Natural birth rate of mosquitoes 

6hm$ Im 

Transfer rate of humans from susceptible to infected compartment 

ah ~ N h 

0hm$ Ih 
N h 

Transfer rate of mosquitoes from susceptible to infected compartment 

Ph 

Natural death rate for humans 

s h 

Death rate of humans due to malaria disease. This death rate in addition to the natural death 
rate is applicable to the humans of Infected compartment. 

A 

Transfer rate of humans from infected to treatment compartment 

€ 

Transfer rate of humans from treatment to recovered compartment 

T 

Transfer rate of humans from Infected to recovered compartment 

<p 

Transfer rate of humans from recovered to susceptible compartment 
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0m.h 

Probability of transmission of infection from an infectious mosquito to a susceptible human 

Ohm 

Probability of transmission of infection from an infectious human to a susceptible mosquito 

r 

Fraction of natural birth rate of humans. Newly born humans enter into susceptible and 
protected compartments with the rates of(l — y)\|/andy\|/ respectively. 

a 

Transfer rate of humans from susceptible to protected compartment 

(A) 

Death rate of mosquitoes. Natural deaths and deaths due to chemicals and other reasons are 
all included. 


Vaccination rate on humans 

£ 

Rate at which Mosquito are killed 

0 

Susceptible mosquitoes bite infected humans with this rate. 

Also, infected mosquitoes bites susceptible humans with the same rate 


Table 1 Parameters of the model (1) and their interpretations 


Variables 

Description 

S h (t) 

Number of humans insusceptible compartment at time t 

P(t ) 

Number of humans in protected compartment at time t 

4(t) 

Number of humans in infected compartment at time t 

T(t ) 

Number of humans in treatment compartment at time t 

R(t) 

Number of humans in recovered compartment at time t 

Sm(t) 

Number of mosquitoes in susceptible compartment at time t 

/ m (0 

Number of mosquitoes in infected compartment at time t 

N h { t) 

Total human population at time t 

Nm(t) 

Total mosquito population at time t 


Table 2 Variables of the model (1) and their interpretations 


1.3. Mathematical formulation of SPITR model 

Applying the assumptions, definitions of compartmental variables and parameters described in tables 1 and 2, the system 
of non-linear differential equations which describe the dynamics of malaria transmission with controlling measures are 
formulated and presented in this section. 

"dS Jdt = (1 - y) y/- 9hm ^ ' m S h - n h S h - gS h + cpR 
d? / dt = yy/+ gS h + (1 - (p)R - \i h P 

dh/dt = ej ^S h - rift - lift - + 5ft)Ift 

< dT/dt = A\ h -eT - g h T (1) 

dR /dt = eT + xl h — (1 — cp)R — cpR — g h R 

dSjdt = p- * h S m ~ uS m 

d\Jdt = 6 -^S m -co\ m 

Proposed model 

We assumed that susceptible and infected mosquitoes are killed at same rate by humans and the inclusion of a vaccination 
parameter gives the modified SPITR model below. Again applying the assumptions, definitions of compartmental variables 
and parameters described in tables 1 and 2, the system of non-linear differential equations which describe the dynamics of 
malaria transmission with controlling measures are formulated and presented in this section. 


(pR 



(1 -<p)R 

Figure 2 Flow chart of SPITR malaria modified model 
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r dS Jdt = (1 -fiy- - fi h S h - gS h + <pR - 

dP / dt = yy/+ gS h + (1 - (p)R — [i h P 
d\ h /dt = ej ^S h - rl h - A\ h - (p h + 8 h )\ h 

< dT/dt = A\ h — eT — \i h T (2) 

dR/dt = eT + xl h - (1 - (p)R — cpR — p h R 

dS m /dt = p - - ^ lh S m - coS m - aS m 
dl m /dt = e -^S m -a)l m -eI m 


The initial conditions of the system of equations (1) are given by 5^(0) = S h0 ,P(0) = P 0 , 4(0) = 4o>T(0) = T 0 , R(0) = 
R 0 ^m(°) = s m 0 and / m (0) = 7 m0 .Also ; we have used in equation (1) that a h = ( 6 mh 0 l m /N h ) and a m = ( 6 hm 0 \ h /N h ). 
The term a h denotes the rate at which the susceptible humans become infected by infectious female mosquitoes. 
Similarly, the term a m denotes the rate at which the susceptible mosquitoes become infected by infectious humans. The 
rate of infection propagated to susceptible humans by infected mosquitoes is dependent on the total number of humans. 
Similarly, the rate of infection propagated to susceptible mosquitoes by infected humans is dependent on the total number 
of humans. 


2. Analysis of the model 

We now analyze the SPITR model in order to show the two controlling methods considered here have substantial impact 
on controlling the spread of malaria disease. In fact, the disease will be completely eradicated if the controlling methods 
are implemented effectively. The two controlling mechanisms proposed here have such a big potential. We consider now 
the solutions of the system of non-linear differential equation (1). We understand that the interpretations of these 
solutions must be biologically meaningful. Hence it is easy to identify that the feasible region of system (1) is M+. The 
seven dimensional solution space shows that all the solutions are positive. Hence, the feasible region containing all the 
solutions of the system of equations (1) is given by the set/2 = {(, S h ,P , I h ,T,R,S m , I m ) E M+}. Here the quantities 
S h ,P , I h ,T ,R ,S m , I m are all non - negatives. Further the total human and mosquito populations are represented by N h 
and N m they have the upper asymptotic values (y//fi h ) and (p/(a) + £)) respectively. Therefore, the region /2 is positively 
invariant i.e. solutions remain positive for all the temporal values. Thus, the model (1) is biologically meaningful and 
mathematical well-posed or well present in the domain /2. 

On summing up all the individual equations from (la) to (le) of the system (1), it is straight forward to get ( dN h /dt ) = 
( y/— g h N h — <44). Here the notation N h = (S h + P + I h + T + R) represents the total human population contained in all 
the five compartments. We consider the solution of the system of equations (1) when the term 8 h \ h vanishes. In case if the 
death rate of humans due to malaria disease is considered to be free, i.e., 8 h = Othen we obtain (dN h /dt) = (y/ — 

The solution of this differential equation is found to be N h (t) = (y//Ph) + [4o — (WPh)] e _M/lt showing that N h (t) -> 
y//li h diSt -> oo. The termA4 0 denotes the initial total human population. It can be interpreted that the total human 
population grows and asymptotically converges to a positive quantity given by (y//ti h ) under the condition that humans 
do not die due to malaria infection. TYmsy//g h is an upper bound of the total human populationA4(t ) i.e.TV^ (°°) < 
y// 1 u h . Whenever the initial human population starts off low below ( y//g h ) then it grows over time and finally reaches the 
upper asymptotic value (^/p^). Similarly, whenever the initial human population starts off high above (y///i h ) then it 
decays overtime and finally reaches the lower asymptotic value (^7^) [19]. 

Similarly on summing up all the individual equations from (If) to (lg) of the system (1), it is straight forward to get 
dN m /dt = p — (a) + c) N m . Here the notation N m = (S m + I m ) represents the total mosquito population contained in all 
the two compartments. The solution of this differential equation is found to be 2V m (t) = (p/o)) + [2V m0 — 
(p/(o) + c))] e -a)t showing that N m (t) -> (p/(ca + c))as t -> oo. The term N m0 denotes the initial total mosquito 
population. It can be interpreted that the total mosquito population grows and asymptotically converges to a positive 
quantity given by (p/(ca + c)).Thus (p/ca) is an upper bound of the total mosquito population^(t ) i.e. N m ( oo) < 
(p/(a) + c)). Whenever the initial mosquito population starts off low below (p/(ca + c))then it grows over time and finally 
reaches the upper asymptotic value(p/(ca + £)). Similarly, whenever the initial mosquito population starts off high above 
(p/(ca + £))then it decays over time and finally reaches the lower asymptotic value(p/(ca + £)). 


2.1. Disease free equilibrium point 

Disease free equilibrium points are steady state solutions when there is no malaria in the human population and there is 
no plasmodium parasite in the mosquito population. That is, absence of malaria causing infections occurs in both 
populations at the disease free equilibrium point. The disease free equilibrium point is denoted by E 0 = 
(S^,I^,P*,T* f R*, 54, Im)- The equilibrium point is obtained on setting the right-hand side of the non-linear system (1) to 
zero. Thus, at the equilibrium point the quantities satisfy the condition = P* = T* = R* = 1^ = 0,5^ = ( y //+ V F)) 
and 54 = (p / ((u) + £))).A1 so, Overhead star represents the values of the functions at the disease free equilibrium point. 
The disease free equilibrium point represents E 0 the disease free situation in which there is no malaria infection either in 
the society or in the environment. Therefore, the diseases free equilibrium point is given by 


E 0 = 


f ~r) ipQ¥+g+Yv h ) 

\n h +g+'v’ Vh^h+g+'V)’ 


0,0,0 ,p/(a) + 


O.o) (3) 
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2.2. Basic Reproduction Number R 0 

Reproduction number, denoted by R 0 , is the threshold or a level for many epidemiological models. It determines whether a 
disease can attack the population or not. The threshold quantity /^indicates the number of new infected individuals is 
produced by one infected individual. When R 0 < 0 each infected individual propagates the infection and produces on 
average less than one new infected individual so that the disease is expected to die out completely over time. On the other 
hand if R 0 > l,each individual produces more than one new infected individual so we would expect the disease to spread 
more and grow in the population. This means that the value of threshold quantity R 0 in order to eradicate the disease 
must be reduced by less than one. We determine the expression for R 0 using the next generation matrix approach [18] 

asR a = — ———. Further, it can be verified that the disease free equilibrium point E a given by (4] is 

0 vKco+£) 2 fa h +g+'¥)(T+A.+ [ i h +S h ') M H 05 y 1 3 

locally asymptotically stable if R 0 < land unstable ifR 0 > 1. 

3. Existence of Backward Bifurcation 

We intend to determine the stability of the endemic equilibrium and to carry out the possibility of the existence of 
backward bifurcation due to existence of multiple equilibrium and reinfection. As a disease attacks it reduces the number 
of susceptible individuals in the population, which tends to reduce its reproductive rate. For a backward bifurcation to 
occur, this means that when R 0 < 0 the endemic equilibrium point can exist as well as when R 0 > 1. We would expect the 
disease to be able to attack at R 0 = 1 in the case of a backward bifurcation with the properties of unstable equilibrium 
bifurcating from the disease-free equilibrium when R 0 < 1, giving rise to multiple stable states. But not in the case of a 
forward bifurcation, in which in the absence of a low-level unstable equilibrium when R 0 < 1 and a stable equilibrium 
bifurcating from the disease-free equilibrium when R 0 > 1 arise naturally when the disease does not attack when R 0 = 
1. A simple criterion for a backward bifurcation, then, is one in which the disease can attack when R 0 = 1. This implies that 
the disease-free equilibrium may not be globally asymptotically stable even if R 0 < 1. 

4. Numerical simulation 

In this section we consider the simulation study of the system of differential equations given in (1). As stated earlier these 
equations describe the dynamics of human and mosquito populations of the malaria model that includes intervention 
strategies. The simulation study is performed using ode45 solver of MATLAB software. The Runge-Kutta fourth-order 
method based on a variable step-size is used for the purpose. The parametric values have been collected from the 
literature and used here. Those were not available were not obtained from literatures published by researchers in malaria 
endemic countries which have similar environmental conditions, we present the numerical analysis of the model(l). Those 
were not available were not obtained from literatures published by researchers in malaria endemic countries which have 
similar environmental conditions, we present the numerical analysis of the model(l) and (60). The initial conditions used 
were S m = 47186, P(0) = 600,4(0) = 47473,7(0) = 1300,R(0) = 47470, S m(0) = 17500,7 m (0) = 26,250. The time- 
axises in all the phase portraits.. We simulate the basic malaria model in the absence of any intervention and the malaria 
model without intervention strategies, and find out the effects of varying each intervention parameter [1],[20],[35] and 
[34]. 

4.1. Estimation of Parameters 

Table.3 We estimate that it will take 3 times a day for 7 days to recovery from malaria infection through 
Chemotherapy and the i ncubation period of malaria in humans was considered [26],[28 ], [20],[2], [22] and [35]. 


Symbol 

Values 

Source 

Y 

0.111 

Mwamtobe, 2009 [20] 

V 

0.000027 

Calculated 

P 

0.0655 

Fekadu,et.al, 2015 [35] 

a h 

0.07143 

Estimated [2 23] and [34 35] 

«m 

0.0909 

Estimated [2 231 and [34 351 

P 

0.0000548 

Calculated 

S h 

0.00000071 

WHO malaria report,2014 [21 and [111 

E 

0.00722 

Gumel,2008 [211 and [221 

g 

0.78 

Mwamtobe, 2010 [22] 

X 

0.01 

Miranda, 2009 [241 

(A) 

1/25 

Blayneh, 2009 [231 

T 

1/7 

FMOH,2004 [251 and [101 

<p 

1/121 

Estimated [2, 23 1 and [341 

@mh 

0.42 

Niger, 2008 [201 

Ohm 

0.0655 

Estimated [35 341 

8 

0.0034 

WHO, 2020 [101 and [121 

¥ 

0.0909 

WHO, 2020 [101 and [121 

0 

0,40 

Chitnis, 2005 [201 and [251 


We also compare the protected and treated classes in order to observe the impact of the interventions on the spread 
control of malaria. Here we consider two intervention strategies viz. INT and IRS, as controlling measures of malaria. 
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Time(days) 

Figure 3 Dynamics of population sizes of human compartments R 0 = 0.8741 

Figure 3 shows the comparative population sizes of five human compartments. Here we have considered two intervention 
strategies viz. INTs and IRS to protect the human individuals of the protected compartment. The situation in the figure 1 
corresponds to the reproductive ratio R 0 = 0.8741. That is, if the intervention strategies are implemented to that order so 
as to reduce R 0 to a value equivalent to 0.8741 the situation shown in this figure occurs. As a result (i) the population size 
of the protected compartment increases, (ii) the population sizes of susceptible, treated and recovered compartments 
decrease and (hi) the population size of the infected class rises slightly for some time and then falls. Further, the 
comparative study promises that the filing of protected compartment leads to emptying the remaining compartments and 
nullifying the malaria disease among the population. The control measures ITNs and IRS reduce the availability of hosts 
and kill mosquitoes. The spread of malaria disease is reduced as the mosquitoes that are attempting to feed themselves 
from human hosts is reduced. There is an increase in the infection when the protection is low and decrease in the infection 
when the protection is high. That is, the infection and protection are inversely proportional to each other. The two 
preventive measures recommended here to protect humans from mosquitoes are effective and economical and hence can 
be implemented. The netting also acts as a protective access to against bites, making it an ideal prevention method in low 
level areas where mosquitoes are found more in size. 



Figure 4 Comparative population dynamics of protected human and infected mosquito compartments when R 0 = 

0.8741 
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Figure 4 shows the comparative population dynamics of protected human and infected mosquito compartments. The 
population sizes of protected human and infected mosquito compartments are proportional to each other. If the protected 
population size of humans increased using intervention techniques then the infected mosquitoes size is decreased and vice 
versa. 

4.2. Simulation of Protected and Treated Human Populations 



Figure. 5 shows the comparative diagram of the protected human individuals and treated human individuals with 

time 

The Figure 5 contains the diagram of the protected and treated individuals. It has been shown that at the beginning stage 
of the treatment has more effective compared to prevention. But with time, the prevention measure plays a bigger part in 
reducing the spread of malaria disease. 



Figure.6 Shows the comparative diagram of the protected and treated human population with time as we vary 

rates X and g 
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5. Conclusions 

We presented and analyzed a basic deterministic 
mathematical model to understand better how to control 
the spread of malaria. Then protected and treated classes 
were added to the basic SIR model and extended it to 
SPITR mode. The protected class is due to the introduction 
of the two intervention strategies ITNs and IRS. The 
potential impact of protection and treatment strategies on 
the transmission of malaria disease are considered. The 
present model includes the protection features those are 
effective to control the transmission of malaria disease 
and used worldwide. 

The numerical analysis of the model suggested that 
effective control or eradication of malaria can be achieved 
by the combination of protection and treatment measures. 
We have that three intervention strategies are combined; 
there is a greater reduction in the number of infected 
individuals. The prevention strategies played a greater 
role in reducing the number of infected individuals by 
lowering the contact rate between the mosquito and 
human populations, for instance, through the use of ITNs. 
On the other hand both prevention strategies led to the 
reduction of the mosquito population hence lowering the 
infected mosquito population. Effective treatment 
consolidated the prevention strategies and also treating 
and exercising vaccination taking part in the pilot of a new 
tool with the potential to save tens of thousands of lives in 
Africa. This paper provides useful tools for assessing the 
effectiveness analysis of a combination of the intervention 
strategies and analyzing the potential impact of 
prevention with treatment. 
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